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ABSTRACT 

The study of the distribution of baryonic matter within dark halos enriches our understanding of 
galaxy formation. We show the radial dependence of stellar baryon fraction curves derived for 21 
lensing galaxies from the CfA- Arizona Space Telescope LEns Survey by means of stellar population 
synthesis and pixel-based mass reconstruction. The sample covers a stellar mass range of ~ 
2 X 10^ — 3 X 1O^^M0 (solar masses) which corresponds to a total enclosed mass range of Ml — 
7 X 10^ — 3 X 10^'^ Mq on radial scales from 0.25i?„ to 5i?„ (effective radii). By examining the Ms 
and Ml dependence on radial distance to the centre of each galaxy we find that there are pairs of 
lenses on small to intermediate mass scales which approach at large radii the same values for their 
enclosed total mass but exhibit very different stellar masses and stellar baryon fractions. This peculiar 
behaviour subsides for the most massive lensing galaxies. All the baryon fraction profiles show that 
the dark matter halo overtakes the stellar content between 1.5 and 2.5i?„. At 3i?e most of the stellar 
component is enclosed. We find evidence for a stellar baryon fraction steadily declining over the full 
mass range. Furthermore, we shed light on the Fundamental Plane puzzle by showing that the slope 
of the Ml(< R)-to-Ms{< R) relation approaches the mass-to-light relation of recent Fundamental 
Plane studies at large radii. We also introduce novel concentration indices c — R9Q/R5Q for stellar 
and total mass profiles (i.e., the ratio of radii enclosing 90% and 50% of the stellar or total mass). 
We show that the value c — 2.6 originally determined by light profiles which separates early- type 
galaxies from late-type galaxies also holds for stellar mass. In particular, less massive dark matter 
halos turn out to be influenced by the distribution of stellar matter on resolved scales below 10 kpc. 
The ongoing study of resolved baryon fraction profiles will make it possible to evaluate the validity of 
star formation models as well as adiabatic contraction prescriptions commonly used in simulations. 
Subject headings: gravitational lensing - galaxies: elliptical and lenticular, cD - galaxies: evolution - 
galaxies: halos - galaxies: stellar content - dark matter. 



1. INTRODUCTION 

The physics driving the evolution from the collapse of 
gas and dark matter halos to the formation of galaxies 
remains one of the open questions in astrophysics. In 
general, star formation efficiency — viewed as the stel- 
lar to total mass fraction within the virial radius of a 
halo — is highest for galaxies similar to the Milky Way, 
with an efficiency decreasing towards higher and lower 
masses (jMpster et al. 2010). The lower escape veloci- 
ties in less massive g alaxies allow t he gas to be ejecte d 
by stellar feedback (ILarsonl [1971 : IDekel fc Silk! [l98fih . 
Supernova-induced winds are energetic enough to sig- 
nificantly impe de galaxy formation at baryonic masses 
below 10"Mo ([Brooks et all [20071 ). Such feedback reg- 
ulates the star formation efficienc y, which is responsibl e 
for the mass-metallicity relation ( Tremonti et al.|[2004[ ). 
For more massive galaxies, an Active Galactic Nucleus 
(AGN) i s believed to account for the decreasing efficiency 
(see e.g. iDi Matteo et al.ll2005l ). This feedback mecha- 
nism explains the exponential cut-off in the luminosity 
function, either by the thermal coupling of AGN out- 
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flows with gas fe.g. iTabor fc Binnevlll993[ : ICroton et al.l 
120061 : iBower et al.ll2006h or bv mechanic al feedback that 
prevents gas cooling (jSiiacki et al.ll2007l ). 

Observational estimates of the stellar baryon fraction 
are thus an essential piece of the puzzle and provide 
important constraints on simulations, especially at 
the sub-grid level that describes the baryon physics of 
galaxy formation. They also help to understand the 
nature of the scaling relations, such as the Fundamental 
Plane and its projections. Currently, most of the studies 
that resolve the central regions of galaxies on scales 
below 10 kpc are based on dynamica l models applied 
to the kinematics of sta rs (see e.g. iCappellari et al.l 
120061 : iCoccato et al.l l2009t ) . Similarly, lensing studies 
on galaxy scales are usually based on a parametric 
decompo sition of the stellar and dark ma.tter c omponent 
(see e.g. lAuger et all 120 lOt iTrott et al.l 120101) . with its 
inhere nt d egeneracies. Over larger scales. iGuo et al.l 
(j2010( ) and iMoster et all (|2010t ) match the stellar mass 
function of SDSS galaxies with the distribution of dark 
matter halos from numerical simulations to find stellar 
baryon fractions /& 3 — 4% — significantly lower 
than the cosra ological fraction — fli,/il„i — 0.17 
punklev et al.l [2001 ) — with a maximum for galaxies 
with halo masses around IO^^Mq. However, this ap- 
proach is only valid for masses enclosed within the virial 
radius, and cannot resolve the radial dependence, which 
offers valuable information about how baryons build 
galaxies. For instance, the velocity dispersion analysis 
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of iLintott et al.l (|2006[ ) on a sample of SDSS early- type 
galaxies gives within the effective radius a low baryon 
fraction (/(, ~ 8%) which is lower than the cosmological 
value, but twice as large as determined within the 
virial radius, illustrating the importance of a resolved 
estimate of the baryon fraction within galaxy halos. 
Galaxy formation models combining the evolution of 
the dark matter and gaseous components along with 
a set of sub-grid p rescriptions for star formation and 
feedback (see e.g. 'Ka uffmann et al.l 119931 : iCole et all 
119941 : iCroton et al. 2006i) are only indirect methods with 
considerable uncertainties. Indeed, robust observational 
estimates of the baryon fraction on galaxy scales are 
needed to properly constrain the recipes included in 
these models. 

Gravitational lensing opens a door to smaller scales 
over which baryonic processes are important. For 
instance, one can explore concentrations and baryon 
fractions giving goo d evidence of adiab a tic co ntraction, 
as done e.g. by iJiang k Kochanekl ()2007D . They 
analyze the relation between stellar baryon fraction and 
concentration in adiabat ic and non-adiabatic models. 
iMandelbaum et all (|2006[ ) present a galaxy-galaxy weak 
lensing analysis of a large sample of early and late-type 
galaxies. They obtain stellar surface masses depending 
on radius with a resolution down to 10 kpc. However, 
their approach is also based on a halo-model to describe 
the relation of galaxies and dark matter. 

The use of mass models and model-based prescriptions 
introduces hard-to-quantify deviations from real mass 
distributions, especially over the scales of a few Re that 
we want to investigate. Assuming a mass model for a 
lensing system excludes mass distributions which are not 
accessible in the parameter-space of the model and in- 
troduces the problem of model non-uniqueness. To avoid 
this, free-form methods are ne cessary. In this paper , 
we use the PixeLens method of (|Saha fc Williamal2004l : 
[Colcs 2008) to reconstruct the surface mass density of 
a sample of lens galaxies. For the stellar component — 
which represents the vast majority of the baryons in the 
inner regions of early-type galaxies — the photometric 
data is used to constrain a larg e volume of stellar 
popu lation synthesis (SPS) models (jFerreras et al.|[2005l 
120081 ). The combination of both lensing and stellar mass 
in a pixel-based manner allows for a two-dimensional 
mapping of the baryon fraction. Choosing a sample of 
moderate redshift lenses enables us to determine the 
lensing profile out to a few Rg. The GASTLES sampl^l 
fulfills this requirement. We present in this paper an 
analysis of the enclosed stellar and total mass content in 
a sample of 21 lensing galaxies out to a radial distance 
of 1.5 — 2 times the Einstein radius, i.e. up to several 

In section [2] we discuss briefly the lens sample with re- 
spect to environment, lens morphology and photometric 
properties. By means of three lens systems, arguably 
rather extreme, we illustrate the subtleties of photomet- 
ric modeling and the authenticity of lenses. The lat- 
ter point refers to unlensed double quasars which mimic 

^ http:/ /www. cfa.harvard.edu/castles 



a lens system with a doubly imaged quasar. We will 
show how a real lens can be distinguished from a spu- 
rious system in our analysis. We test the reliablility of 
our photometry-based results by comparing inferred stel- 
lar surface ma ss densitie s with equivalent results from 
iFerreras et all (2009. ) and lShen et all (pOOl . 

Section U) presents the results of this study regarding 
the radial dependence of enclosed stellar versus total lens 
mass. We continue in section |5| with a closer examination 
of the stellar and total mass concentration and define a 
simple model to study the energetic evolution of early 
type galaxies. The conclusion and discussion section [6] 
summarizes our findings and puts them into the context 
of recent work on galaxy formation. 

2. SAMPLE PROPERTIES 

In the following we compare lens samples in general 
with respect to their environment. We also explain how 
the environment affects the lens systems and continue 
with a detailed one-by-one study of the lenses used in 
this analysis according to their photometric and morpho- 
logical properties. To get a clear view on lens galaxies 
whose baryonic content we want to determine, we need 
to correct the data for light originating from the quasar 
images, by means of PSF subtraction and masking. 

The selection criteria for our lens sample are as follows. 
For the lens mass reconstruction we must have the red- 
shifts of the source and the lens as well as accurate image 
positions. The stellar population synthesis analysis de- 
mands a sufficient separation between lens and quasar 
images in order to extract uncontaminated photometric 
estimates from the lens. Furthermore, NIR imaging must 
be available. 

Constraining the SPS models using photometry in sev- 
eral bands is desirable, although we note that our refer- 
ence H-band is the F160W fiher of HST/NIGMOS. For 
the redshift of most of the lenses, this band maps a rest- 
frame region that does not change significantly for the 
colours found in these galaxies. We discuss in this section 
the available multiband data and respective PSFs used 
for the modeling of the surface brightness distribution. 
Finally we discuss outliers and special cases for compar- 
ison. All information regarding lens galaxy properties, 
their environment and photometry is given in Table [T] 
and Tabled 

2.1. The Environment 

To describe a lens with respect to its environment, 
one has to keep in mind that the lens shear required 
by (parametric and non-parametric) lens models can be 
due to physically interacting galaxies or to line-of-sight 
objects. Regarding the former, one could estimate how 
the environment of the lens galaxy evolved in its recent 
past, whereas any line-of-sight objects are naturally un- 
related to the local region of the lens. Nevertheless these 
two sources for shear are hard to distinguish. If located 
in a group or cluster environment. X-ray measurements 
are exp ected to gi ve reliable constraints on the DM 
content (jBuote fc T sai 1995) and thus a hint about the 
direction and strength of the shear. Only a few lens 
environme nts have been studied so far for CASTLE S 
lenses (e.g. lFassnacht et al.|[200llMomcheva et al.ll2006D . 
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The environ ment for a sample o f 70 SLACS lenses has 
been studied bv lTreu et all (|2009D : they find 17±5% are 
in overdense regions. For our sample of 21 CASTLES 
lenses we find that 7 galaxies are located in groups and 
3 in clusters. Four galaxies have one close galaxy or pos- 
sible companion with which they may interact gravita- 
tionally. For the remaining 7, no large shear contribution 
is required and no close galaxies have been found. Thus 
we find that ~ 50% of our galaxies lie in overdense re- 
gions. The lower fraction found bv iTreu et al.l (|200 9f) is 
likely due to the smaller redshift range of SLACS (up to 
z 0.5) and the property of the SDSS selection func- 
tion to pick lenses whose Einstein radius is about the 
fibre-radius of the SDSS spectrograph (3 arcsec). 

2.2. The Sample 

In the following we briefly discuss the lens sample. In 
addition to the previous paragraph we provide infor- 
mation about the environment which in fact influences 
both the mass model and the light profile and yields 
important insights into the evolution of early-type 
galaxies. After giving the full name of the lenses we use 
their abbreviations only. The method used to model the 
lenses is explained in detail in section 13.31 Figure [15] 
shows the appearance of the free-form mass models. 

For 9 lensing systems all three bands were used to 
constrain a large number of SPS models, which conse- 
quently sets constraints on the colour-to-mass relation. 
Another 8 lenses could be analysed in H and I band. 
The remaining 4 lenses had suitable data in H band only. 

First, we describe the 9 lenses with suitable data in all 
three wavebands. 

The four-image lens system ("quad") B0712+472 is 
one of the few lenses for which a TinyTim PSF was 
sufficient to remove quasar images in H band. In V 
and I band the quasar images could all be masked out. 
Lens models found in previous studies require signifi- 
cant external shear, which can be attributed to 9 or 
mo re g alaxies in a for eground group at z ^ 0.3 found 
by iFassnacht fc LubinI (2002) . Their study also shows 
one other galaxy at the redshift of B0712 at ^ 100" from 
the lens. 

The quads B1422+231, B2045+265, Q0047-2808, 
Q2237+030 undergo the following treatment. In both 
I and V bands, TinyTim provided a suitable PSF. In 
the H band an isolated star taken from the same or a 
contemporaneous NICMOS image was used for convo- 
lution and point-source fitting if needed. B1422 is in a 
poor group with 5 nearby galaxies mostly south east of 
the lens that cause a significant shear (fMonichc va et al.l 
I2006t iHogg fc Blan dford 19 9j). The group is v isible in 
X-rays at 0.5 - 2 keV (Momcheva et al.l IBM ). In re- 
cent work by I Wong et all (|20lH ) 12 new members were 
found to be part of the group. B2045 as found by 
IFassnacht et al.l ()1999() might be infiuenced by a group 
of galaxies west of the lens. A shear in this direction 
is also required by the lens model. The lens might also 
be affected by a close dwarf ga laxy causing anomalous 
fiux ratios (jMcKean et al.ll2007D . Q0047 is a lens with 
only a smal l shea r required by lens models. However, 
iWong et al.l (|2011f) find evidence for a galaxy group with 
9 members. In the case of the Einstein Cross Q2237 the 



bulge of a spiral galaxy is responsible for the lensing. The 
system shows only a mild external shear due to the disk 
of the spiral galaxy. 

For spiral galaxies the contribution of dust to the pho- 
tometry is usually more significant than for early-type 
galaxies (the latter morphological type constitute the 
majority of our lensing galaxies). However, we note that 
in the case of Q2237, the redshift of the lens is very 
low, which implies that our reference photometric band 
(H) maps a similar wavelength range in the rest frame, 
where d ust attenuation i s less severe. From the esti- 
mates of lEigenbrod et all (pOOl on VLT/FORSl spec- 
tra of Q2237, we infer a contribution from dust in the 
H band photometry of Q2237 at the level of 0.05 mag 
(jFerreras et al.ll2010 V 

For lenses BRI0952-0115, Q0142-100 and 
PG1115+080 extensive use of the iteration method 
described in section 13.1! was made if the quasar images 
could not be masked out. The environments of the 
doubly im aged quasars BRI09 5 2 and Q0142 have been 
studie d bv iLehar et al.l (I2000D . iMomcheva et al.l ()20Q6f ) 
and Eigcn brod et al. n2007f ) and found to have no 
dominant impact on the total shear beyond a cosmo- 
logical (large-scale structure) contribution 7lss which 
is additionally confirmed by lens models. BRI0952 was 
previously thought to reside in a regi on loosely bound 
to a poor group with 5 members (|Momcheva et al] 
12006( 1: a later study found it is at hig her redshift and 
thus not connected with the group (Eigcnbro d et all 
|2007| ). For Q0142 there is not much known about 
the close-in group environment, although there are 
some galaxies near the line-of-sight , whose redshifts 
are mostly unknown. iSurdei et al.l ()1987D speculate 
that a galaxy about 10" away from the lens may be a 
group member. The env ironment of the q u ad PG 1115 is 
thoroughly analyzed by IMomcheva et al.l (|2006f ) . They 
find 13 galaxies in a local gro up with eloiigated group 
emission m X-rays according to lGrant et"all (|2004D . The 
brightest 4 members of the group are located on an axis 
with a position angle of +60° (measured North through 
East) of the lens mass which accounts well for the shear 
required in our lens model. 

The two-image lens HS0818+1227 requires spe- 
cial treatment as we use an isolated PSF of image 
B1030-f071 to fit the quasar image in the H band. In 
the I and V bands, the quasar images are used for 
fitting. Iteration as it is used for enhancing PSFs of 
other lenses does not provide better model fits for the 
lens because of the large separation between images 
and lens. The image separation is 2.56". Hence the 
reduction process is further simplified by masking. Since 
its discovery by Hagen & Reimers (2000) no further 
insights into the environmental properties o f the l ens are 
available. Nevertheless Hagen fc Reime rs! ()2000( ) found 
a galaxy 5" north of the lens which appears to have the 
same redshift oi z = 0.39, which explains the external 
shear required by our lens model. A chain of galaxies 
at a distance of 10" north-east could also be associated 
with the lens galaxy. 

Next, we describe the 8 lenses with suitable data in 
two wavebands. 

For the quad B1608+656 and the doubles HE1104- 
1805 and HE2149-2745 H and I band data could be 
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used t o isolate the lens galaxy. According to iMore et alj 
(|2009f) MG2016+112 exhibits quadruply imaged fea- 
tures of the quasar jet which can be distinguished only 
in the radio band. We take account of the rather complex 
structure in the lensing part of our analysis. For all their 
H band images, a sufficiently isolated star with fitted 
background extracted from the image of MG0414-I-0534 
was used to remove the quasar images with an acceptable 
goodness of fit. 

B1608 resides in the middle of a galaxy grou p with 
8 oth er group members according to Fassnach t et al.l 
(|2006( ). The photometry shows an object close to the 
main galaxy, which constitutes a second lens galaxy. This 
is confirmed by the reconstructed mass map (Fig. [151) 
it predicts a conspicuously elongated mass distribution 
towards NE. Images also show a prominent dust lane be- 
tween the two galaxies. However, we analyzed the impact 
of dust reddening on our results, as shown in Appendix 
[B]for B1608 and B1600. The uncertainty due to dust on 
\og{Ms) is in both cases not larger than 0.3 dex. 

MG2016 is known to be a giant elliptical galaxy in a 
cluster with 69 probable photometrically selected mem- 
bers of many different galaxy types (Toft et al. 2003). 
Among them is a significant fraction of merging cluster 
galaxies, which is direct evidence for a h ierarchical for- 
mation history ()van Dokkum et"alll2000[ ). Most of the 
neighbouring objects within 30" lie on an east-west axis 
and thus explain the major shear direction. HE1104 fea- 
tures the second highest image separation of 3.19" and 
a distinct lens galaxy (the median separation is ~ 1.5"). 
Furthermore the lens appears to be near the bright im- 
age which is rather unusual and implies the presence of 
a gro up or cluster enhancing the separation (jLehar et al.l 
12000'). Parametric as well as free-form mass models also 
suggest that an external shear is mandatory to repro - 
duce the image configuration (e.g. IWisotzki et aLlll998D . 
The lens galaxy is unaffected by quasar light allowing for 
a good fit. However, the photom etric redshifts of a few 
neighbouring galaxies described in lFaure et al.l ()2004l ) in- 
dicate that such cluster galaxies are probable companions 
of the lensed quasar rather than of the lens. The dou- 
ble HE2M9_jinght_be a member of a cluster as inferred 
by iLopez et aLl ()1998[ ) by a large number of red non- 
stellar objects in R-band images of the field around the 
lens. Considering recent es timates of the lens redshift 
from lEigenbrod et al.l (120071) Czi.ns = 0.6 03) and the en- 
vironment survey from IMomcheva et al.l (2006) HE2149 
could be in a group with 3 neighbouring objects. The 
morphology of the lens shows no sign of strong external 
shear. 

The doubly imaged quasar SBS1520+530 is treated 
like the previous doubles but with a star from the same 
image file in preference to other PSFs. This lens is a 
member of a. galaxy group with at least 4 other members 
as stated in lAuger et al.l (|2008f ). 

For the two quads MG0414+0534, 
RXJ0911+0551 and the double Q0957+561 we 
obtain good residual maps by means of the iteration 
method. MG0414 at z = 0.960 is the second most 
distant lens of our sample. Judging by its luminos- 
ity and colour, the lens is likely to be a passivel y 
evolving early-type galax y (|Tonrv fc Kochanea 119990 . 
iSchechter fc Moora (|1993f) find an object close to image 
B visible only in I-band, which might contribute to 



the lensing effect. Our reconstructed mass map also 
shows increased surface density at the position of 
the obje ct. RXJ0911 is loc ated on the outskirts of a 
cluster (j Morgan et al.l I2001|) . Chandra observations 
of the cluster suggest a complex non-spherical cluster 
mass distribution a t a temperat u re of r oughly 2.3 keV. 
Q0957, found by iWalsh et al.1 (|1979f ). is special in 
several ways. First there is a doubly imaged galaxy 
component in addition to the famous double quasar used 
to calculate the projected mass map. Secondly the lens 
is a cD galaxy located in the centre of a cluster. The 
nearest cluster member lies within 10" East of the lens 
galaxy. However a simple external shear is insufficient 
to describe the effect of the environment on the image 
positions. Breaking the degeneracy between the shape 
of the galaxy and the cluster she ar takes advantage 
of arc features (Keeton et al. 2000) and X-ray data as 
attempted bv Chartas et al.. (1998.) . 

Finally, we describe all the lenses with suitable data in 
only one waveband. 

B1030+071, B1152+200 and B1600+434 are 
treated similarly with regard to the fitting routine, i.e. 
the isolated outermost quasar image was used for sub- 
traction and convolution. The three doubles have compa- 
rable angular image separations and average velocity dis- 
persions as well as intermediate luminosities. Observed 
substructures in B1030 i ndicate the presence of an in- 
teracting galaxy system (jJackson et al.l [21300^ although 
firm statements ab out the environment cannot be made 
(jLehar et al.ll2000t ). However, shear is not strongly re- 
quired by our mass model. For B1152 there is no infor- 
mation about the composition of the environment. Judg- 
ing by the morphology of the image-source system no 
strong shear is expected. B1600 is located in a denser 
group with at least 6 late-type galaxies (jAuger et al.l 
[2007D which cause significant shear. The absence of X- 
ray emission is suggestive of a not relaxed group, a con- 
clusion strengthened by the elongated morphology of the 
group. Furthermore the lens galaxy appears to be al- 
most edge-on and exhibits a prominent dust lane. As 
remarked above, dust reddening changes the population 
synthesis input and leads to underestimated stellar con- 
tent, but even for the extreme cases in our sample the 
effect of dust on inferred stellar masses cannot be larger 
than 0.3 dex (see Appendix IB)) . 

For the doubly imaged quasar LBQS1009-0252 the 
star in the H-band image of MG0414 is u sed again as 
a PSF with sufficient quality for the fit. iLehar et al.l 
(|2000D locate the lens galaxy close to quasar image B. 
They find that a dominant shear contribution of the host 
galaxy of a n earby quasar (4.6" northwest of the lens — 
unrelated to the lensed quasar) is consistent with the 
derived major axis of the lens when modeled by a Sin- 
gular Isother mal Ellipsoid. Using a singular isothermal 
sphere model | Claeskens et al.l (|2001[ ) determine a smaller 
shear. iFaure et al.l ()2004f ) state that there is no signifi- 
cant galaxy overdensity in the field. This is in agreement 
with the free-form lens models of this study, which do not 
require external shear for this lens. 

2.3. Outliers And Special Cases 

We now briefly describe three special cases, 
B0218+357, B1933+503 and RXJ0921+4529. 
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With the first two we want to demonstrate the impact 
properties hke small image separations and interfering 
luminous structures can have on the goodness of the 
SPS. The third lens shows how spurious lenses, i.e. 
galaxies with nearby quasars which are not lensed 
images of the same background object, behave in this 
analysis. All three lenses are excluded from our analysis. 

For B0218 as for 10 other systems in our sample a star 
was used to fit the quasar images in the H-band. Since 
B0218 is the system with the smallest image separation 
(0.33") known, it is extremely difficult to separate 
the lens galaxy from the images of the background 
quasar. The system is an extreme case in several 
aspects and a good example for showing the impact 
of degeneracies between the magnitudes of overlapping 
objects. B0218 unlike any other lens in the sample 
did not yield reasonable Sersic profile parameters as 
the wings of the quasar PSFs overlap with the lens. 
For an unconstrained fit the combined light from the 
quasar images and lens galaxy results most likely in an 
overestimated magnitude of the PSFs. However after 
attempting to fit the lens system only by PSFs, a Sersic 
profile is needed to achieve an acceptable residual map. 
Even though one cannot obtain zero residuals by fitting 
only two point sources, there are several combinations 
of Sersic profile magnitudes and two PSF magnitudes 
that result in the same total surface brightness profile. 
Bearing this in mind, we use the fitting parameters with 
the best which also yields an acceptable residual 
map, to carry out the SPS. The projected lens mass 
map shows that external p otentials ind u ce a shear in 
B0218 that was studied in iLehar et all (|2000f ). They 
find 13 possibly perturbing galaxies inside a radius of 
20" located roughly along the axis which connects the 
two quasar ima ges. It should be ni entioned that B0218 
is according to iLehar et al.l ()2000[ ) a late-type galaxy 
which causes the SPS to predict a different mass content. 



B1933, discovered by iSvkes et al.l (|1998D . has 10 
distinct images formed from a three-component source, 
promising an exceptionally well-constrained mass pro- 
file. A star in the same H-band image was used for 
convolution. There is as yet no study of the environment 
of the lens but according to the mass reconstruction, 
no strong shear is necessary to explain the morphology. 
The resolved features of the lensed background object 
cannot be fitted by PSF but are taken out of the fitting 
routine by using circular masks with a 5 pixel radius, 
a size chosen to cover features distinguishable from 
background and still show enough of the lens galaxy to 
allow for a reasonable fit. The trade-off between light 
contamination due to minimal masking and information 
loss due to aggressive masking is in any case problematic. 
In the case of B1933 almost the whole inner region is 
surrounded by masked regions causing the fit parameters 
i?e and n to diverge. Setting a constraint on the Sersic 
index {n < 4) is necessary. Despite all attempts at 
modeling this lens, it remained a persistent outlier, and 
hence is removed from the analysis. 

The double RXJ0921, has the highest angular image 
separation (6.93") compa red to any oth er lens in the 
sample. According tO iMuiioz et al.l (|2001l ) it is probably 



a member of an X-ray cluste r. From model fits of the 
host galaxy iPeng et ai.l (|2006[ ) conclude that RXJ0921 is 
a binary quas a r rath er than a gravitational lens. Also 
iPopovic et al.l ()2010l ) find quite different spectral prop- 
erties in the spectra of the two components. For now 
we assume the system is a lens. Since even the smaller 
lens-image distance is above 3" and the quasar images 
are isolated, we obtain a high-quality fit by taking the 
quasar image as a PSF for both overall convolution and 
quasar subtraction. No constraints are necessary. There 
are 16 objects within 20" from the lens galaxy. Only 
for three of them a redshift close to that of the lens 
could be determined. The mass model however does 
not require an external shear. In contrast to all other 
lenses, RXJ0921 (when treated as a lens) turns out to 
exhibit an unusually low stellar-mass fraction and an al- 
most constant Ml{< R) profile. The peculiar properties 
of RXJ0921 can be taken as further evidence against the 
lens hypothesis as suspected in aforementioned studies. 

3. ANALYSIS TECHNIQUE 

In this section we describe necessary steps to obtain 
baryon-fraction profiles. As explained before we need to 
free the lens galaxies in given multiband data from in- 
terfering light, originating from quasar images, to obtain 
best-fitting surface brightness profiles. This reduction 
step is shown in the following paragraph. The output 
is used to constrain SPS models and to estimate pixel- 
by-pixel the stellar mass via the colour-to-mass relation, 
as described in section [3.21 Combining stellar mass esti- 
mates and pixel-based mass reconstruction (section 13. 3p 
yields baryon-fraction profiles for the given set of lens 
galaxies. 

3.1. Preparing Photometry 

The problems arising during this procedure can be as- 
signed to one of the following categories: (a) finding an 
appropriate PSF for convolution and point-source reduc- 
tion, or (b) removing perturbing light sources that neg- 
atively affect the fitting procedure. The latter includes 
masking of image regions as well as fitting of additional 
light profiles to structures which clearly do not belong to 
the lens. Since we use photometric data in different fil- 
ters (V,I and H bands), one of the following PSF-picking 
procedures has to be suitably chosen for each band. 

1 . Find an isolated star from the same or a contempo- 
raneous (as nearly as possible) image and extract it 
and a sufficiently large surrounding region not con- 
taminated by light from the lens system or other 
sources; hereafter referred to as the star-picking 
method. 

2. Select the outermost image of a lens system and 
use it for the quasar image fitting. While the 
lens galaxy and the other quasar images are fitted 
with GALFIT (jPeng et al.ll2002[) . the residual im- 
age, showing only the previously chosen outermost 
image without any contaminating light, can then 
be taken as a qualitatively refined PSF. This step 
can be repeated until we reach the desired level of 
enhancement. This procedure will be referred to 
as the iteration method. In some cases, when the 
picked quasar image was sufficiently isolated, iter- 
ation brought no further improvement. 
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3. Using a sy nthetic PSF generated by TinyTim 
(|Kristl 119931) yields better results for WFPC2 im- 
ages rather than for NICMOS data (presumably 
due to the higher stability of WFPC2 PSFs). 

Methods (1) to (3) were combined with masking of 
quasar images or any other luminous structure interfer- 
ing with the fit. To prevent the fit from diverging, in 
some cases further constraints are necessary. 

Details of the constraints applied to each lens can be 
found in Table [2] In the following, we give an overview of 
the types of constraints we have applied to the sample: 

• fixing the sky background for the already reduced 
images to a value we determined with SExtrac- 
tor, since estimating the background is essen- 
tial to extract a m eaningful profile of the lens 
(|Haussler et al.ll2007l ). (1 ow signal-to-noise objects 
are thus neglected, increasing the goodness of the 
fit for the generally bright lens galaxy), 

• fixing the surface brightness profiles to previously 
determined (x, y) positions, 

• constraining i?„ and/or the Sersic index n, since 
both parameters are degenerate, being basically in- 
versely related, i.e. constraining i?„ to a low value 
causes n to diverge and vice versa, 

• constraining the position angle and the axis ratio to 
a physically appropriate range of values, and finally 

• restricting the range of magnitudes of the point 
sources, e.g., constraining image A to be at least 
0.5 magnitudes brighter than image B. 

Except for the first two constraints, which are neces- 
sary for only a few lensing systems, we try to keep the 
number of degrees of freedom of the fit as high as possi- 
ble and fix the parameters only if there is no alternative. 
In cases where the best-fit PA in one band was found 
to differ significantly from the PA in another band, it 
was necessary to constrain that parameter. The same 
applies to the axis ratio b/a. If a highly eccentric ellipse 
is fitted to an actual round lens galaxy due to interfer- 
ing PSF wings, the parameter space must be constrained 
to exclude less likely b/a values. Since both the lensing 
galaxies and quasar images are in some cases too bright 
to be distinguishable in H-band, but too faint in V-band, 
we use the shape parameters PA or b/a from I-band as 
a proxy for the fits in other bands. This was done for 
B1422, BRI0952, HS0818 and Q2237. This approach is 
legitimized by the fact that different stellar populations 
visible in different bands do not change their relative po- 
sitions and orientations considerably. All other fitting 
parameters apart from the boxiness, which is set to zero 
throughout the process are free. T he boxiness as wel l as 
all other parameters are defined in lPeng et al.l (|2002| ). 

To minimize and test the stability of the fit, the 
fitting procedure was repeated with slight changes to the 
initial parameters. But is not the only criterion to 
assess the quality of a fit. We focus on the goodness of the 
fit in regions most important to our analysis, i.e. in the 
central region of the lens galaxy. The ratio of the fitted 
image and the original image yields a percentage map 
of the lens systems showing pixel-by-pixel the quality 




Fig. 1. — Upper left panel: original H-band NICMOS image of 
the lens system B2045. North is left and East is down. Best fits 
for the lens galaxy could be obtained by masking out images A, B, 
C, the stellar objects south of the lens and a "blob" West of the 
lens galaxy. Upper right panel: model for the lens galaxy. Lower 
left panel: Residual image. Lower right panel: Residual divided by 
original image. 

of the model, as illustrated in Fig. [T] We note that 
in some cases, better fits were rejected in favour of 
the flatness of the residual in central regions < li?„ of 
the lens. This occurs in particular for the few lenses for 
which the short distance to a lensed quasar image makes 
constraints mandatory. 

For lenses with photometric data in more than one 
band the I or V-band parameters for R^, n, b/a and 
PA were taken as a prior to the H-band parameters if 
necessary. 

3.2. Estimating stellar masses 

Our stellar mass estimates are based on a pixel-based 
comparison of the best fits to the surface brightness of the 
lenses with stellar mass-to-light ratios (T) determined 
by population synthesis models constrained by the avail- 
able photometry. Even though GALFIT does a para- 
metric search to get the best fit, for this analysis we 
are just interested in the 2D distribution that minimises 
the residuals, regardless of the parameters themselves, 
i.e. we are less sensitive to the inherent degeneracies as- 
sociated with parametric fits. For each lens we ran a 
grid of 32 X 32 X 32 models, where the star formation 
history is described by a decaying exponential, defined 
by three free parameters — the quantities in parenthe- 
ses denote the range explored for each one: formation 
epoch (defined as a redshift 2 < zpoR < 10); exponen- 
tial timescale (—1 < log(T/Gyr) < 1); and mctallicity 
-1 < [m/H] < +0.3). Mo dels from l Bruzual fc CharloJ 
200l are used, assuming a lChabrieii (|2003( ) initial mass 
function (IMF). For each choice of the three parame- 
ters, a composite population is obtained, transformed to 
the redshift of the lens, and folded with the passband 
response of the HST V (F555W, WFPC2), I (F814W, 
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WFPC2) and H (F160W, NICMOS) filters to compare 
with the observed colours and to extract a mass-to-light 
ratio in the observer-frame H band. The colours are 
corrected for Galact ic extinction using the dust maps of 
ISchlegel et aP (1998). 

We extract the stellar mass densities from the H band 
image (NICMOS F160W). Whenever model fits of the 
lens were available for I or V, the colours where used on a 
pixel by pixel basis to constrain Th. Otherwise, we used 
integrated colours within an elliptical aperture defined 
by the half-light radius of the H band image (see Ta- 
ble [1} . In general, broadband photometry alone cannot 
be used to constrain the ages and metallicities of the lens 
galaxy. However, the stellar masses, when estimated via 
"red" M/L ratios, ar e less sensit ive to the age-metallicity 
degeneracy (see e.g. lFerreras_et al. 2008). For compar- 
ison, we provide the total stellar mass-to-light ratio in 
the rest- frame V-band Ms/Ly in Table [TJ Colours and 
m agnitudes are in agr eement with comparable quantities 
in lRusin et all (|2003l) . The F160W band corresponds to 
a rest- frame wavelength between 0.8 and 1.2fim (except 
for Q2237, which roughly samples rest- frame H-band). 
Hence, for the sample considered here, the mass-to-light 
ratios are not affected by the presence of young stars, an 
issue that becomes importan t when dealing with optical 
or NUV indicators (see e.g. iRogers et al.l l2010l) . From 
the modeling of the old stellar populations that these sys- 
tems feature (except for lens Q2237 — which is a bulge 
— the other lenses are early- type galaxies), an uncer- 
tainty of AT < 0.15 dex is expected (Gallazzi fc Belli 
I2009D . Dust reddening, as we explain in Appendix jB] 
leads to underestimated stellar mass. However, since no 
starburst galaxy is among our lenses we can safely as- 
sume that the effect of dust on Mg does not exceed 20%. 
The number of lenses which exhibit dusty features (e.g. 
B1600) is, nevertheless, small. The most significant sys- 
tematic error relates to the choice of the Initial Mass 
Function, especially the low-mass end, which does not 
contribute to the light, but can contribute very signif- 
icantly to the total mass conten t. However, frequently 
used c hoices of the IMF such as iMiller fc Sca lo (1979), 
IScaiol (Il986h . IKrouna et all (I1993D or IChabrier (2003i) 
have similar distributions at the low ma ss end. It i s only 
the traditional single-power law of the iSalpeteii (|1955f ) 
IMF that gives different stellar mass predictions. Pre- 
vious detailed w ork on the kinematic s of nearby early- 
tvpe galaxies f i Cappellari et al.l 120061 ) or strong lenses 
(jFerreras et al.ll2008L l20ioi r shows that the low-mass end 
of the Salpeter IMF is ruled out as it predicts stellar 
mass surface densities higher than the dynamical or lens- 
ing estimates. However, even using a Salpeter IMF does 
not strongly affect our results, as they mainly focus on 
the scaling of t he reg ions where dark matter dominates. 
iFerreras et al.l (|2010l ) illustrate differences between five 
different population synthesis models based on different 
prescriptions and/or stellar libraries. The predicted stel- 
lar masses - measured in the H band - agree to within 
10%, (at fixed IMF) especially given the ages of these 
lenses. 

To compare our lensing (early-type) galaxies with a 
typical field sample, we show in Fig. [5] the equivalent to 
the Kormendy relation (this time defined with respect 
to the surface stellar mass density at li?„). We show as 
open dots the sample of ACS/GOODS early-type galax- 
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Fig. 2. — Stellar Surface Mass density at a radius of 1-Ro versus ef- 
fective radius. Early- type galaxies with stellar mass above 10^° Mq 
(open circles) selected from the Hubble Space Telescope/ Advanced 
Camera for Surveys images of the Great Observatories Origins 
Deep Survey (GOODS) are shown together with our lens sample 
(filled squares). The dashed line denotes the stellar mass to size re- 
lation from SE)SS which accounts for early- type galaxies at ^ ~ 0.1 
llShen et al.|[2003.') . 

ies from (jFerreras et al.|[2009() . and our lensing galaxies 
as squares with error bars. One can see that 14 out of 18 
lensing galaxies are located inside a lu-band around the 
best fit of the GOOD S sample. We also provide the SDSS 
relation from lShen et ai.l (|20Q3f ) as a local {z 0.1) refer- 
ence. The obvious preference of the lens sample to be at 
larger effective radius and smaller surface mass density 
is due to a selection bias, which is a combination of the 
lensing bias and additional requirements, such as a suffi- 
cient distinguishability from surrounding quasar images. 

Furthermore, as we show later in this paper, the slope 
of the Fundamental Plane relation ^ L can be re- 
covered from our data. Thus we consider the lens sample 
representative for early-type galaxies in general. 

3.3. Reconstructing the total-mass profiles 

For each lens, the projected total-mass distribution is 
reconstructed on a circular field made up of 750 square 
tiles or pixels, each pixel consisting of a uniform non- 
negative mass distribution with a mass density of a few 
times the critical density. We provide mass reconstruc- 
tion maps of the sample in Fig. |T5] The pixellated mass 
distribution must reproduce the observations, in the fol- 
lowing ways. 

1. Multiple-image systems with the observed posi- 
tions must arise as solutions of the lens equation. 
The images are considered to be unresolved; for ex- 
tended images, the peak of their surface brightness 
distribution is located and considered as an unre- 
solved image. Note that although the mass distri- 
bution is discontinuous at the pixel boundaries, the 
lens equation is continuous. 

2. For lenses with measured time delays, the model is 
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required to reproduce them. Concordance values 
of Ho,nM,^A (i.e. 72 km s"^ Mpc^\0.3,0.7) are 
assumed. 

In addition, the mass distribution must satisfy the fol- 
lowing prior conditions. 

1. The local density gradient must point no more than 
45° away from the centre of brightness. Since the 
central regions of galaxies are expected to be dom- 
inated by stars, it seems safe to assume that the 
mass and light peaks coincide. 

2. The circular average (around the centre) of the 
projected density, falls off as R^^^^ or faster. 
The three-dimensional mass profiles of galaxies are 
thought to be invariably steeper than r~^'^, so 
again this appears to be a safe prior assumption. 

3. No pixel is allowed to be more than twice the mean 
of its neighbours, except for the central pixel, which 
can be arbitrarily high to mimic central density 
cusps, 

4. Unless the lens shows signs of asymmetry, the mass 
distribution is required to be symmetric under a 
180° rotation around the centre. 

In practice, there are infinitely many mass models 
that satisfy all the above conditions, becau se solutions 
of the lens equation are highly non- uniq ue (jFalco et al.l 
[19851: [Sahal[2000t ILiesenborgs et al.l [20081 ). Accordingly, 
for each lens we generate an ensemble of 300 models, by 
a random-walk technique in model space. The random 
walk implicitly defines a prior measure in model space, 
and it turns out that an equivalent statement of the prior 
measure is that it m ust be invariant under rescalings 
of units (|Colesl [20081 ). A very useful property of the 
above constraints (observational or prior) is that they 
can all be formulated as linear equations or inequalities. 
This means any weighted mean of ensemble members is 
also an admissible model. Hence, we can conveniently 
use the ensemble mean to represent a typical model. 
However, the main results later in this paper use the full 
model ensembles. 

The lens models do not attempt to subtract off lens- 
ing mass outside the galaxy. Such mass could come from 
the environment, 14 of the lenses being in dense envi- 
ronments (see column env in Tab. [T|), or it could be in 
an interloper along the line of sight. However, given that 
the model-ensemble technique yields conservatively large 
error-bars on the mass maps, we expect that external 
lensing mass is unlikely to be larger than the estimated 
uncertainties. 

In the following section, we will consider the circularly- 
averaged enclosed mass profile M{< R), see Figures [3| 
and m The outermost radius to which the mass pro- 
files are reconstructed is fixed to two times the lensing 
radius -Ri„„b, which is defined as the radial position of 
the outermost lensed image with respect to the center 
of the lens. We choose 2i?,„„g as a trade-off between un- 
certainty and common radial range for the sample. The 
range of enclosed-mass profiles in the ensemble, which is 
interpreted as the uncertainty, has a characteristic but- 
terfly shape. That is to say, M(< R) is well constrained 



in the image region, but becomes more uncertain farther 
in or out. Note that the butterfly shape is less promi- 
nent or even distorted for less symmetric lensed image 
configuration. The steep limit of the butterfiy shape is 
expected to be roughly M{< R) ~ R^-^, resulting from 
the minimal steepness of in the projected density. 

The shallow limit of the butterfly shape is given by the 
steepest model in the ensemble. 

4. RADIAL DEPENDENCE OF STELLAR VERSUS TOTAL 

MASS 

To compare the radial dependence of stellar and to- 
tal mass, it is interesting to consider pairs of lenses with 
matching Ml{< R) or with matching Ms{< R). To il- 
lustrate, in Fig. 131 we show three pairs of galaxies with 
the following properties (see also Table [T|): 

1. small mass, matching Ml{< R) profiles, differing 
Ms{< R), 

2. intermediate mass, matching Ml{< R) profiles, 
differing Ms{< R), 

3. high mass, differing Ml{< R) profiles, matching 
M,(< R). 

The radial scale is R/Rsin where i?Ein has been esti- 
mated from the pixelated mass maps. Error bars are 68% 
confidence from the population-synthesis models used for 
Ms{< R) values. For Ml{< R) we use error bars cor- 
responding to 90% of the Ml{< R) r ange of the model 
ensemble, as described in section 15751 Note that the er- 
rors attached to Ms{< R) and Ml{< R) are correlated. 

The matched pairs are, of course, only rough matches. 
Also, the i?Ein values are not the same for the matched 
pairs of galaxies. The angular radial scale is proportional 
not only to the enclosed mass but also to {d^s I {dLds))^'^ 
(corresponding to {d^dLs / ds)'^'^ for the physical Ein- 
stein radius), where the d's are the angular diameter 
distances between observer and lens (L), observer and 
source (S) and lens and source (LS). Latter distance ratio 
must be approximately equal to identify matching pro- 
files. To enable comparison between scales we include 
Rsin/R, in Table [H With these caveats, we point out 
some interesting features. 

Consider first the two low-mass lenses PG1115 and 
Q0047 (bottom panels. Fig. [3]). While the total mass 
within 2i?i„„5, is very similar, the stellar mass of PG1115 
rises only to 50% that of Q0047. The same qualitative 
behavior is seen if these two galaxies are compared using 
rather than R/R^^^ as the radial scale. Neverthe- 
less, these two low-mass lenses have fs = Ms/Ml ^ 0.17 
the range of high baryon fractions. Lenses within this 
range can consequently be referred to as high /s lenses. 

Comparing the two intermediate mass lenses B1030 
and MG0414 in the middle panels of Fig. [31 we find 
that their cumulative total mass curves are very similar. 
However, the stellar mass of B1030 is just ~ 30% that 
of MG0414, independent of the radius. If we consider 
the stellar radial scale, we find that MG0414 has ^ 4 
times the stellar mass of B1030. Their baryon fractions 
approach values from /s 0.05 (B1030) to fs ~ 0.17 
(MG0414) at the outermost radius to which we have es- 
timates. In the intermediate mass range of our sam- 
ple (roughly 5 x lO^M© to 15 x 10"Mq) MG0414 has 
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Fig. 3. — Cumulative stellar mass and lens mass profiles against 
projected radius (in units of the Einstein radius) for three compara- 
ble pairs of lenses. PG1115 and Q0047 (bottom row) are low mass, 
B1030 and MG0414 (middle) are intermediate mass, and HE1104 
and RXJ0911 (top) are high mass. The gray vertical line marks 
2i?e- To assist comparison, each right-hand panel has the profiles 
from the corresponding left-hand panel duplicated with thin dotted 
lines. 

one of largest and B1030 the lowest stellar-mass fraction. 
It should be noted that the opposite behaviour, namely 
matching stellar profiles on both i?ie„, and scale with 
very different total mass is also possible. An example for 
the latter would be a comparison between B1030 (mid- 
dle row) and PG1115 (bottom row) of Fig. [3] with equal 
Ms{R/ Re) but total mass profiles differing by a factor of 
4.5 at 3i?e- 

The two high mass lenses RXJ0911 and HE1104 have 
a total stellar mass of ~ 2 x lO^^M©. For comparison, 
the total mass profiles differ slightly for radii < l.hR^ and 
> 2. Si?,,. But it should be noted that RXJ0911 is located 
at the centre of a galaxy cluster which might lead to a lens 
mass estimate slightly larger than the actual virial mass 
of the lens galaxy. At 6.5i?„, i.e., 2i?,„„, HE1104 has 
12% less total mass than RXJ0911. At 3i?„, i.e. i?,„„, 
the difference is still 6%. In terms of stellar-mass fraction 
HE1104 exhibits small values of fs w 0.07 and RXJ0911 
of fs w 0.06. In the high mass regime (> 15 x IO^^Mq) 
the range of possible stellar-mass fractions appears to 
be small compared to low and intermediate masses, and 
always close to 0.05. Those two lenses are thus represen- 
tative for low fs lenses. 



Comparing lens profiles on R^ scales intrinsic to the 
luminous part of the galaxy, one can find many lenses 
with similar stellar mass profiles, which is not surpris- 
ing. After all the enclosed mass values Ms{< 2Ri^„^) 
cover with ^ 10^" to 2 x 10^^ a relatively small 
range in contrast to a total mass range Ml{< 2_R,,,„J 
of ~ 2 X 10^° to 2 X 10^^. However, pairs of lenses 
with matching Ms{< R) and Ml{< R) profiles over the 
whole radial range as shown in Fig. |3] are rare. Most 
lenses with matching Ms{< R) profiles exhibit quite dif- 
ferent Ml{< R) profiles. The above lenses HE1104 and 
RXJ0911 are - apart from their data points > 6R^ - 
matching pairs within uncertainties on i?„ scale as they 
are on i?,„„s scale, a consequence of i?io„s/i?c being equal 
for both objects. 

In Fig. U we present the left column lenses of Figl3] 
now on baryonic scales, two of them have new counter- 
parts with similar Ms{< R) and Ml{< R). As before we 
present low to high mass galaxies from the bottom up. 

For the two low mass lenses PG1115 and B0712, we 
find that at their outermost common radius ^ 2.5R^ 
their baryon fraction is - 0.08. B1030 and BRI1009 also 
match well within their error bars although the mean 
stellar mass profile of B1030 is consistently below the 
one of BRI1009. The error region of its lens mass profile 
shows quite large error bars and thus make it easy to 
match. The baryon fraction at 2.7 R^ is approximately 
fs — 0.08. If we compare lenses along the vertical direc- 
tion of Fig. |3]B0712 and BRI1009 are representative for 
most lenses on low to high mass scales, that is, similar 
Ms{< R), dissimilar Ml{< R) and baryon fractions. 

In summary we find on both baryonic scale R^ and 
lensing scale Ri^ns- 

• many pairs with the same enclosed total (lens) 
mass, but with different enclosed stellar mass, 

• a small number of pairs (decreasing with increasing 
Ml) with the same enclosed total (lens) and stellar 
mass. 

We can already conjecture an anti-correlation between 
enclosed lens mass and stellar-mass fraction, which will 
be studied in detail later on. 

However, one should keep in mind that our result could 
be influenced by the lens environment and its history. 
See also Table [1] column ^Env^ and section [2] for infor- 
mation on the local lens environment. The phenomenon 
of same Ml but different Mg becomes less prominent for 
larger total lens masses, on both i?,e„g and R^ scale. Nev- 
ertheless, global trends and interdependencies might be 
revealed by analysing the whole set of lenses, which is 
done below. 

Using our sample of 21 lensing objects we consider the 
following relations to highlight the interdependencies in 
the {Ml,Ms,R) parameter space: 

1. the enclosed total mass Ml{< R) as a function of 
enclosed stellar mass Ms{< R) at a fixed radius R, 

2. the stellar- mass fraction as a function of radial dis- 
tance, fsiR) = Ms{< R)/Ml{< R), 

3. the stellar-mass fraction as a function of the total 
mass AIl, 
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Fig. 4. — As in Fig.[3l but with the effective radius as reference 
scale, shown for PG1115 and B0712 (bottom), B1030 and BRI1009 
(middle), HE1104 and RXJ0911 (top). Here the grey vertical line 
marks -Ricns. 

4. the stellar-mass fraction as a function of redshift. 

Fig. [5] shows the first relation for a range of radial 
positions from 0.25R^ to 5R^, parametrized by the di- 
mensionless quantity x = R/R„. For reference, we list 
in Table [T] the enclosed stellar and lens mass within 2i?„ 
with error bars. 

The universal baryon fractio n according to WMA P5, 
fb = fls/^M = 0.17 ±0.02 (IHinshaw et al.ll200l . is 
included as two dashed lines for the upper and lower 
bounds. The solid line denotes a stellar-mass fraction of 
one, i.e. the total mass content consists of 100% stellar 
mass. Note that the data points refer to baryonic matter 
in stars and do not account for other baryonic content 
like gas and dust. The gas contents of our lens sample - 
mostly early type galaxies - is expected to be small. But 
for the Einstein Cross (Q2237), which is the bulge of a 
spiral galaxy, and B1600, which is likely to be a late- type 
galaxy viewed edge-on, one can indeed expect deviations 
from the obtained Mg values. 

The galaxy B1608 shows an unreasonably high stellar- 
mass fraction for radii < .75i?e (e.g. left panel in top 
row of Fig. [SJ. To take proper account of the light dis- 
tribution we fit both the brightest galaxy and its merging 
companion with Sersic profiles, but we only use the infor- 
mation of the light profile of the brightest galaxy for the 
computation of stellar mass. The enclosed mass values 



are thus taken with respect to the center of the brightest 
galaxy. As a consequence of the degeneracy between the 
two Sersic profiles the central region of the light profile 
is modeled rather poorly. This causes an overestimate of 
the stellar content (< 15%) in a region where the neigh- 
bouring galaxy, which is also responsible for light deflec- 
tion, interferes with the fit. The pixels with highest total 
mass and highest stellar content do not match for B1608. 
This also causes larger deviations in the region < li?„. 

The late-type galaxies Q2237 and B1600 might be sub- 
ject to dust reddening. In general the impact of redden- 
ing on high redshift lenses is stronger due to the bluer 
populations observed in H-band and the higher absorp- 
tion of dust at smaller wavelength. However, on the basis 
of the analysis shown in Appendix |B1 we do not expect 
departures of more than 20% towards higher Ms- This 
will shift i?1600 closer to the bulk of lenses in Fig. [5] 

The prominent fs curve of B1422 — starting at twice 
the value of most other lens galaxies — might also be 
caused by light contamination. This time it originates 
from the innermost quasar image which lies just 0.25" 
away from the galaxy centre, an extreme among the 21 
lenses studied in this paper. 

In the online material of this paper we provide a movie 
version of Fig. [5] to help visualize the trend of the stellar- 
baryon fraction with increasing radius. See also Ap- 
pendix 1^ 

The lens galaxies reveal the following properties, which 
are qualitatively assessable already from Fig. [5l but will 
be explained in detail later on: 



1. Most lenses populate a band of 0.1 < fs 
within 5i?,,. 



< 0.4 



2. The slope of the enclosed Ml-Io-Ms relation of Fig. 
[S] within the shown radial range becomes gradually 
steeper for larger enclosed radii (an effect quanti- 
fied in the following paragraph). 

3. Between 2 (1.5) and 2.5i?e (2i?c) for most lenses 
with total mass below (above) 4 x IO^^Mq the dark 
matter halos overtake the stellar content, that is 
they move primarily toward increasing total mass. 
The turning point thus depends on the halo mass. 
The dark matter halos of more massive galaxies 
start to dominate the matter balance at larger radii 
(in units of i?„) than those of less massive galaxies. 

Note that by "overtake" we refer to the radius where 
dAl^/dR K. dMs/dR rather than to the radius where 
the total stellar mass contributes 50% of the total mass. 
As a consequence of limited resolution this radius can 
only be given with larger uncertainties (^ 0.5-Re)- 

Point 2 can also be illustrated by plotting the slope 
77 determined from Mg oc so that it represents light 
as a function of mass. We find that rj asymptotically 
approaches 0.75, as one can see in Fig. [6l which is in 
agreement with previous studies of the fundamental 
plane within error b ars (e.g., Guzman et al.l 119931 : 
iJorgensen et al.l 119961 : iLeied l2009l )^ A bootstrapping 



method for a large and a reduced sample is used to 
determine the Mg-to-Mi relation and its standard errors 
respectively. Both runs are done with 10^ realizations. 
The 19-lens sample contains all the lenses except for 
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Fig. 5. — The panels show the enclosed Ml against enclosed AIs plane for a number of apertures, defined by the radial distance xR^ to 
the centre of the lens galaxy, labelled by 'x' in the top-left corner of each panel. We cover a radial distance from 0.25 X ij^ to 5 X i?o from 
upper left to lower right panels in conveniently chosen steps. Grey circles highlight a subset of 8 lenses which are probed out to 5-Ro. The 
solid line denotes the equality of total and stellar mass, whereas the dashed lines represent the upper and lower limit of the global baryon 
fraction IHinshaw et al.ll2009l ). 



the outlier B1608 and the late-type galaxy Q2237. 
Farther out in radius, the number of lenses with profiles 
extending to a particular radius decreases. Because 
of that, Fig. ini also shows the number of lenses used 
for each fit. As a consequence of changing sample size 
discontinuities appear between 2.25 and 3.5i?„ and at 
4.5i?o- The most extreme ones are caused by B0712 
(2.5i?,) and B1030 (4.5i?J falling out of the sample. The 
behaviour of the error bars in a bootstrap fit depends 
on the size of the drawn sample subset. To get more 
meaningful error bars we fixed the size of the sample 
subset to be 50% of the available number of lenses at 



each radius. The small sample instead comprises all 8 
lenses being probed out to 5i?e, which are highlighted 
in Fig. E] by grey filled circles. From M, oc M^"^^^^-^^ 
at 0.25i?e the reciprocal slope rj{R) declines as 1/R 
and ends up at 5i?„ with the relation Ms oc 
We expect only small deviations from this slope for 
larger radii since we run out of stars, and additional 
mass from the dark matter halo shifts the distribution 
upwards, whereas possible baryonic contributions from 
gas shift the whole population farther to the right of 
Fig. [5] Additionally, for the 19-lens sample a weighted 
best fit for 77 (i?) suggests that the function approaches 
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Fig. 6. — Slope of the Ma-to-M^ relation taken from Fig. [5]plot- 
ted against the distance to the centre in terms of effective radii. The 
median slopes are determined via a bootstrapping fitting method 
with 10'' realizations to compute meaningful standard errors for a 
sample of 19 lenses (filled squares) and a reduced 8 lens sample 
(open squares). The numbers at the filled squares give the num- 
ber of lenses probed out to the respective radius. The dotted line 
represents a weighted best fit of rj{R) ~ l/ij+const. 

asymptotically a constant value of 0.77 ±0.01. Note that 
for the small sample declines rapidly to reach the 

value of 0.76 already at ~ 1.5R„ and thereafter shows 
no significant departure from it. However, the change in 
slope from small to large radii is significant for both the 
19-lens and the 8-lens sample. 

All the stellar-mass fraction curves in the left and 
right hand panel of Fig. [7] turn over to a similar 
stellar- mass fraction between 1.5 and 2.5 R^, a fact 
also reflected by ri{R) in Fig. [6] With increasing 
radius, the stellar-mass fractions of high mass galaxies 
(Mi(< 2i?J > 4 X IO^Mq) tend towards lower values 
in the majority of cases, meaning fs < 0.2. Low mass 
galaxies (Ml(< 2i?e) ^ 4 x lO-^^M©) show a larger range 
of possible stellar-mass fractions at high and low radii, 
which are in a range between 0.1 and 0.35 (see left hand 
panel of Fig. [7]). This is the reason for the large scatter 
of enclosed stellar-to-total enclosed masses at small radii 
in Fig. [SI 

Averaged over the whole lens sample we find that the 
stellar-mass fraction declines with increasing radius from 
its value /s(< li?o) enclosed in IR^ to only ~ 71% of 
fs,max at 2i?„ ~ 55% at 3i?e, 39% at 4i?„ and finally 
~ 33% at 5i?e- Splitting the sample with respect to total 
mass as done before yields a different picture: For lenses 
with MLi< 2i?J < 4 X lO"Af0 79% of the stellar-mass 
fraction at li?„ remains at 2i?,, 63% at 3i?„, 47% at 4i?„ 
and finally 40% at 5R,. For lenses with Ml(< 2R,) > 
4 X 10 ^^Mq 64% of the stellar-mass fraction at IR^ is 
found at 2i?„, 48% at 3i?„, 33% at 4i?„ and finally 27% 
at 5i?e- The uncertainties of stellar-mass fractions at IR. 



for low Ml lenses are only as high as 10%. For larger 
radius and mass the fs errors decline strongly to less than 
1%. From this we can conclude the following. 

1. Low mass galaxies show a shallower decline in 
their enclosed stellar-mass fraction than high mass 
galaxies: either their stellar content is less concen- 
trated than in high mass galaxies or their dark 
matter content is more concentrated. This point 
becomes clearer in Section [5l where we calculate 
concentration indices of stellar and total mass pro- 
files, 

2. The relative stellar-mass fraction of high versus low 
mass galaxies is significantly offset by a constant 
value within 5i?„ from the centre, i.e. 



fsi< R) 



fsi< R) 



/,(< li?J ^Mr.>4EllMe l^J IM^<4E11AU 



for 2i?„ < i? < 5i?e. 



- 0.15 
(1) 



The latter phenomenon becomes more evident when 
plotting the stellar- mass fraction at fixed R/R^ against 
the total mass as in Fig. [8l From left to right the panels 
show the fg-Mi, relation at discrete radii of 0.5, 1.0, 2.5 
and 4.0 R^. It should be emphasized that the solid line 
fit does not imply a physical relation extendable to the 
high or low mass end of the plot. Note that the relation 
has a tendency to steepen gradually towards lower radii 
whereas the scatter increases. Comparing this to recent 
results from iGuo et all ()2010f) where the ratio of total 
enclosed stellar mass and halo mass Mhaio is analyzed 
with an abundance matching method, we find that their 
stellar-mass fraction curve shows a peak at a halo mass 
of around 6 x 10^ ^M© and decreasing fractions towards 
lower and higher halo masses. This is overplotted in the 
last panel of Fig. [8] 

The effective height of the curve /{, is reduced in con- 
trast to our results owing to the fact that there is signif- 
icant dark matter in the halo extending up to the virial 
radius, which is defined as 



GM, 



halo 



V100iJ2(z) 



1/3 



(2) 



i?vir is roughly a hundred times larger than the region 
probed in this study. We list R^^, values deduced 
from Ms of this study given their Afs-to-Mhaio relation 
in Table [U which has additional implications on the 
lens environment, provided that the lens behaves like 
an SDSS-Galaxy plus simulated halo counterpart of 
respective stellar mass. To visualize how the computed 
stellar-mass fractions change between our resolution 
range and the virial radius we multiply a constant factor 
by the stellar-mass fraction curve from IGuo et all ([20lol ) 
and divide its total mass by the same factor (here we 
use 5). The slope of the high mass end of their curve 
agrees with our best fit of ]\/f£0.i6±().04 .^^j^jjjj^ error bars. 
Scaling to lower radii makes the mismatch for lower Ml 
even more prominent. However, we co nclude that down 
to a certain level the Ms-to-Mhaio from lGuo et al.l (|20100 
is scalable. In the 5i?e-to-i?vir-i'ange, the lower-mass 
lensing galaxies need to decrease their stellar-mass 
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Fig. 7. — Left panel: Stellar-mass fraction against radius in effective radii for lenses with lens mass enclosed within 2Ra below 4 X 10^^ 
Mq. Right panel: Similar, but for lenses with A'/i,{< 2R^) above 4 X 10^^ Mq. 



fractions by a larger amount than hig h mass gal a xies i n 
order to match with the results from lGuo et alj ()2010D . 
Note that scahng our lens sample instead towards higher 
total masses and lower stellar-mass fractions yields the 
same result. 

We should point out that th is direct comparison of our 
results with iGuo et all ()201Cl[ ) is imperfect, since Mhaio 
and Ml are differently defined and the spatial distribu- 
tion of dark matter in a region not directly addressed in 
either paper is unknown. On the other hand the steep- 
est part of the total mass profiles is already enclosed 
and the cumulative mass profiles saturate, i.e. the slope 
of the ML-to-fs relation is only slowly changing beyond 
5i?e. These different trends for fs at lower masses could 
be indicative of an underestimated stellar-mass fraction 
towards smaller halo masses or an overestimated bary- 
onic content towards higher halo masses. If the afore- 
mentioned study of the /s-to-Mi dependency is correct, 
then our findings give rise to the question of what makes 
the stellar-mass fraction of low mass galaxies decline less 
strongly within 5i?o than in the range from 5i?e up to the 
virial radius, in contrast with high mass galaxies. Ex- 
pressed in terms of stellar mass content we find a steeper 
decrease of stellar-mass fractions towards larger Ms than 
our results predict. 

The above defined virial radius i?vir becomes smaller 
for lower stellar mass content. Low Mg galaxies reside 
in halos with larger fa than high Ms galaxies, meaning 
the mass in the dark matter halo relative to Ms is even 
larger, i.e. small galaxies have more concentrated dark 
matter halos than larger ones (see also Section [5]) . 

In order to investigate the influence of the 
distance/lensing-bias, we also show the redshift de- 
pendence of the stellar-mass fraction in Fig. [9l The 
ordinate might be subject to several biases. The lensing 
galaxies plus halo must be massive to produce an 
observable signature. The galaxy should not be too 



faint to be seen and has to obey our selection criterion 
of sufficient separation from quasar image. Fig. |9] shows 
that the correlation between stellar-mass fraction and 
redshift becomes more pronounced with larger radius. 
However, the strongly increasing scatter below AR^ blurs 
the correlation and the slope shows no uniform trend. 

5. DIAGNOSTICS OF BARYON COOLING 

We now consider two different measures of the stellar 
and total-mass profiles, with a view to gaining insight on 
the evolution of lens galaxies from formation to observa- 
tion redshift. 

5.1. Concentration index 

Our spatially resolved stellar and total mass maps al- 
low us to study the difference in concentration of the 
baryon and the total mass d istribution. We defin e a 
concentration index (see e.g. iBershadv et al.l |200C1[) as 
c = -R90/-R50, where i?90 and i?50 denote the radii en- 
closing 90 and 50 percent of the mass (either stellar or 
lens mass). For the luminous component, a concentration 
index above 2.6 indicates an early-type galaxy, whereas 
indic es below 2.6 refer to late-type galaxies(see e.g. Fig. 
1 in iFerreras et al.|[2005l) . Previous studies based on the 
surface brightness distribution use the Petrosian radius 
(or a given number of Petrosian radii) to define the to- 
tal brightness. In our case, we redefine c and take the 
respective radii of our cumulative stellar mass and total 
mass profiles instead — 100% corresponding to enclosed 
masses at 2i?,,„^ (except for Q0957 and IIS0818 where it 
is 1.5i?ie„J. In Fig.[TO]we show concentration versus red- 
shift in the left-hand panel with no obvious correlation 
and the frequencies per concentration bin in the right- 
hand panel. Note that defining the concentration values 
using i?Ein instead of i?ie„s will change the concentration 
values slightly, but even for Ri^^s/Rmn — 1-5 we obtain 
changes in the lens mass concentration of less than ^ 30% 
and only for lenses with high concentrations. 
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- The stellar-mass fraction determined at 0.5, 1.0, 2.5 and 4.0 against total mass. The best fits are found for the sample 
lenses with mean stellar-mass fractions above 1 at R^, which is the case only for B1608. 
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Fig. 9. — The stellar-mass fraction determined at 0.5, 1.0, 2.0 and 4.0 Ra against redshift. 



From Fig. [TU] the frequency distribution of cm^ 
peaks between 3.0 and 3.5 which is in agreement with 
most concentration studies of ea rly- typ e galaxies (e.g. 
lYamauchi eFal] 120051: fPeng et al.ll2010f ). That is, even 
with the redefined concentration quantities one can dis- 
tinguish the lens galaxies morphologically. Evidence is 
given by the two late-type galaxies Q2237 and B1600 
which indeed lie below 2.6. For the two merging galaxies 
in the lens B1608, the interfering potential (for cml) or 
light (for Cm, ) , causes the concentration values to be de- 
creased, pushing CA/g down to 2.6. For the same reason 
we obtain rather large error bars on the lens mass. 

One could check in detail now if the strong correla- 
tion between c and Hubble-type is maintained for the 
newly defined cm^ ■ If '^6 define a concentration parame- 
ter by means of the total mass profiles we expect, as our 
findings in Section |4] already suggest, a totally different 
distribution. Most lenses exhibit cm^ values in a narrow 
region between 1.5 and 2. However, neither in cm^ noi' in 
cml can a clear evolutionary trend be found. Figure [11] 
shows that the concentration parameter for stellar mass 
h^-s a rising trend with total lens mass, whereas cml 
clearly declines with lens mass. 

Note that the error bars of cm^ and cml are the stan- 
dard errors of the R90/R50 values of each model in the 
ensemble multiplied by student's t for a 95% confidence 
interval. This was done since the ensemble can be seen 
as being part of a normal population. The horizontal er- 
ror bars are the Ml errors at the outermost radius of the 
reconstructed mass profile. As we can see at low total 



lens masses the distributions of Mg and Ml are almost 
the same, which means that the Ml profile approaches 
the distribution of the baryonic matter. An interaction 
between the baryonic and dark matter distribution seems 
to be a reasonable explanation, since already in Section 
|4] we find that the stellar-mass fractions of less massive 
lenses are larger than for the more massive lenses. 

A possible interaction between baryons and dark mat- 
ter is likely to influence the slope of the total mass dis- 
tribution close to the center of galaxies. If we assume a 
density following a pure power law p(r) ~ the en- 



closed mass becomes M{< R) 



/(3-^). Thus the 



concentration c and the density slope obey the relation 



/3 = 3 



lnO.9/0.5 

Inc 



(3) 



Figure [12] contrasts the relation between /3 and c based 
on a pure power law (solid line) and data for different 
radial extents. The (3 values represent weighted best fits 
to lens mass profiles with standard errors from the fit. If 
the mass distribution does not follow a pure power-law 
(i?90/i?50) might depend strongly on the radial extent 
of the lens {2Ri^„^). Therefore we compare in both pan- 
els of Fig. [12] concentration values inferred from differ- 
ently sized profiles, with a maximal radius of li?ic„a and 
2i?i„„j,. The innermost data point has rather large un- 
certainties and deviates in most cases strongly from the 
trend at larger radius. To demonstrate its impact on the 
relation we contrast fits with (left panel) and without 
(right panel) regard of the innermost point. According 
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Fig. 10. — Concentration index c = R90/R50 versus rcdshift, 
where c is the ratio of the radii enclosing 90% and 50% of the total 
stellar mass (cm^ i open circles) and lens mass {cm^ , filled squares). 
For MG2016 R50 cannot be calculated due to lack of data points at 
small radii. The solid line indicates c = 2.6 separating early-type 
(c > 2.6) from late- t ype galaxies (c < 2.6). Error bars for index c 
can be found in Fig. 1111 

to eq. Owe find that with increasing Ml{< 2R^) - i.e. 
decreasing concentration - the slope /3 gets shallower. 
It is remarkable how extraordinarily well the weighted 
power law fits neglecting the innermost point reproduce 
the simple /3(c) model at low concentrations, but fail to 
do so at large concentrations, where the data lies below 
the pure power law relation. Higher values of /3 corre- 
spond to shallower profiles. Including the innermost 
point always flattens Ml{< R) fits, which explains why 
respective /? values, although less representative for the 
outer part of the profiles, are in better agreement with 
eq. 131 We can conclude that 

1. excluding the core region of the lenses, we ob- 
tain power law indices and concentration parame- 
ters (i?90/i?50) indicative of a pure power-law be- 
haviour for small concentrations. This notion is 
strengthened by only small shifts of {R90 / R50) go- 
ing from 2 to lRi^„^ profiles. 

2. For more concentrated total mass distributions, we 
find evidence for a significant departure from pure 
power law behaviour. This is confirmed by signifi- 
cant shifts of (i?90/i?50) while reducing the extent 
of the lens from 2 to li?ic„s and an increasing /3 
error towards higher concentrations. 



5.2. Energy ratio 

By means of the stellar mass content one can approach 
the subject of galaxy formation from a different view- 
point. The first question is: "is it feasible to determine 
a characteristic quantity which gives us the amount of 
energy lost between the collapse of an initial sphere of 
homogeneously distributed baryons and its later state as 
a lens galaxy?" . One could ask as well for a ratio of the 
radius of the pre-collapse sphere and an observable spa- 
tial quantity, like the effective radius. Even though this 
is a rough estimate, one can gain insight in the evolution 
process of galaxies. 




Fig. 11. — As in Fig. llOl but plotted against the total mass M^^tot 
enclosed in 2i?c. The y-axis error bars represent for cm^ the stan- 
dard errors for mean concentrations of 300 models. For cji/^ the 
error bars correspond to the uncertainties originating from a 10% 
error in the flux per pixel. 
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Fig. 12. — Power law index /3 of the density profile plot- 
ted against the concentration parameter c = {i?90/-R50)M£, for 
the total lens mass distribution. Filled (open) squares indicate 
c = (iJ90/-R50)M£, determined for a profile extending to 2iJion5 
(IRiens). The left panel uses power law fits to the whole radial 
range, whereas the right panel neglects the innermost data point. 
The solid black line shows the pure power law case according to 
eq. m 

At the time of collapse ii, a region decouples from 
the expansion of the surrounding universe. The baryons 
which are assumed to be homogeneously distributed in 
this sphere are for now assumed to make up the whole 
stellar content of the later lens galaxy, neglecting any 
kind of active evolution such as caused by mergers, ram 
pressure, tidal stripping, etc. The radius ri of such a 
sphere at ti is 



1/3 



(1 + zi)- 



(4) 



where VLb = 0.0441±0.0030 is the baryonic energy densit 



in terms of critical density according to iHinshaw et al 
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Fig. 13. — The stellar-mass fraction determined at 1.0, 2.0, 3.0 and 4.0 iJc versus E2/E\. The solid line denotes the best fit of a power-law 
with slope a. The fits exclude Q2237 and MG2016. 
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Fig. 14. — The energy ratio versus enclosed stellar and enclosed total mass at 1.0, 2.0, 3.0 and 4.0 R^. 



(|2009f l and pc ~ US-STMo/kpc^. The average New- 
tonian energy per unit mass at ti consists only of the 
potential energy per unit mass, which is 

GMtot 



El 



(5) 



(6) 



The total mass Mtot is defined here as 

As the collapse goes on the baryons start to fall in to 
build a more tightly bound structure. At the observation 
redshift, i.e. at time i2, we find mostly objects which are 
in virial equilibrium, that is E = — T, where T denotes 
kinetic energy. Thus we can determine the total energy 
per unit mass of the galaxy at ^2 to be 



where 



Eo = -To = -l<Tf^. 



EE VGAfi(< R)/R 



(7) 



(8) 



is an effective velocity dispersion inferred from lensing 
(see Tabled]). It is computed at i? = i?i„„, and assumed 
not to vary drastically with radius. This effective disper- 
sion has been shown to be an appropriate surrogate for 
the observed kinematic velocity dispersion ((Leier .2009^ . 
Thus, 



El 

E2 



2GM,, 



(9) 



using Eq. SJ 



E2 
E'i 



M. 
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(10) 



Thus we get a quantity E1/E2 oc x m'^^^ /Ml or 

— 1/3 

RicnsMs fs- This is reminiscent of the Kormendy re- 
lation, except that it relates to three-dimensional rather 
than projected densities. For definiteness, we assume a 
formation redshift zi = 5, but the value only implies a 
multiplicative constant. Plotting the energy ratio against 
the stellar-mass fraction, we find a strong correlation 
(Fig. [TH) regardless of the enclosure radius. The slope 
changes only marginally, but the scatter decreases with 
increasing radius. 

However, E2/E1 appears to be uncorrelated with the 
enclosed stellar and total mass. For different radii one 
obtains Fig. [T31 The fact that E2/E1 exhibits such a 
tight correlation with but no clear correlation to con- 
tributing masses, can be interpreted as insensitivity of 
the star formation in early-type galaxies to active evolu- 
tion processes over the time span from zi to ziens- 

6. DISCUSSION 

A resolved, model-independent and thus non- 
degenerate (w.r.t. Ms and Ml for fixed fs) estimate of 
stellar versus total mass within galaxy halos is crucial to 
constrain current galaxy formation models and prescrip- 
tions of baryon-dark matter interactions used therein. 
Besides dynamical methods to explore scales below 10 
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kpc the combination of strong gravitational lensing and 
population synthesis used in this paper is most promising 
to give robust estimates of stellar-mass fractions. 

The analysis of the radial dependence of the mass pro- 
files of 21 CASTLES lenses presented in this paper al- 
lows us to draw the following conclusions. The relation 
between basic galaxy properties, i.e. Ml, Ms and R^^ns 
cannot simply be scaled with their mass. The scatter in 
this parameter space turns out to be especially large for 
galaxies of smaller size. The study of Mg versus Ml and 
of the stellar-mass fractions (/s = Ms/Ml) enables us 
to discriminate between lens galaxies below and above 
Ml{< 2i?J = 4 X 1O"M0. The high mass class pop- 
ulates a lower and narrower /s regime (0.05 to 0.2) on 
given scales and runs out earlier of stellar mass (i.e. at 
lower enclosed radius) than low mass lenses. The latter 
exhibit a more inhomogeneous behaviour with a wider 
range in fs (0.1 to 0.5) and respective slopes. 

We conclude that between 1.5 and 2.hR^ dark matter 
halos start to dominate the matter balance depending on 
their total mass. This M^-dependence causes high mass 
galaxies to gain mass primarily in the form of dark mat- 
ter already at lower radii than low mass galaxies. There- 
fore the slope of the mass-to-light relation, which is a 
projection of the fundamental plane — or our equivalent 
representation, M^J oc Mg — becomes shallower with in- 
creasing radius and asymptotically approaches a slope of 
Tj = 0.76 ± 0.07. Thus the FP tilt can be recovered as a 
gradually growing process with radius. Equivalently, the 
stellar-mass fraction shows a strong correlation to the 
total mass. As we contrast /^(Mi) with a c omparable 
curve deduced by abundance matching from iGuo et aLl 
()2010[ ) dissimilarities for low Ml galaxies become more 
evident the smaller the enclosed region gets. This is likely 



to be a result of different halo definitions, physical prop- 
erties and processes, like baryon-dark matter interactions 
and adiabatic contraction which is beyond the scope of 
the aforementioned study. However, the /s-to-Mhaio re- 
lation scaled down to 4i?„ agrees quite well with lenses 
with Ml ~ 1O^^M0, since the biggest part of stellar 
matter is still enclosed. 

Another important result of this study addresses the 
concentration of stellar (cm, ) and total {cml ) mass pro- 
files. The rule-of-thumb delimiter of c = 2.6 which 
separates early- type galaxies (c > 2.6) from late- types 
c < 2.6) holds also for the concentration parameter (cj\/^ ) 
defined by means of stellar mass instead of luminosity. In 
the low mass regime Ml{< 2R,) < 4x lO^^M© both, cm, 
and cml: tend to similar values around 2.6. This means 
that the total mass profile is likely to be infiuenced by the 
distribution of baryonic matter in stars. From IO^^Mq 
upwards, cm, and cml diverge, due to a stronger confine- 
ment of stars in more massive dark matter halos. The 
cml values above 4 x 10^-^ Mq stay around ~ 2 instead. 
Studying the interdependency of density slope and cml 
wc find that the reconstructed lens profiles show devia- 
tions from a pure power-law mass model, which is evi- 
dence for the sensitivity to the radial trends of the dark 
matter distribution. 

The results presented in this paper are critical to ongo- 
ing studies about the reliability of parametric lens models 
and prescriptions used in galaxy-formation models. The 
radially resolved stellar to dark matter fractions should 
thus also serve as benchmarks for future simulations. 
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APPENDIX 
ANIMATED RESULTS 

An animated version of Fig. [5] is provided in the online material of this paper. The movie contains three panels. 
The left panel shows the enclosed Ml against enclosed Mg plane depending on the size of the aperture, defined by the 
radial distance xR^ to the centre of the lens galaxy. The solid black line denotes the equalit y of total an d stellar mass, 
whereas the dashed lines represent the upper and lower limit of the global baryon fraction (jHinshaw et al. 2009). The 
upper right panel shows the lens PG1115 which is highlighted by a red label in the left panel. The lower right panel 
shows stellar baryon fraction versus radius as in Fig. [71 The solid black line denotes the baryon fraction curve of 
PG1115. With each time step of the movie the enclosure radius increases indicated by the factor x in the legend of 
the left panel and the red lines in the two right-hand panels. We cover a radial distance from 0.125 R^ to 5 R^ in 40 
time steps. 



DUST REDDENING 

In this section we study the impact of dust reddening. Fig. shows what happens when dust is included in the 
analysis of two lenses where the contribution of dust could be important: B1600 and B1608. For this exercise we take 
a simpler set of models, but the underlying effect from dust is similar to the more complex case involving exponentially 
decaying star formation histories. We run a set of simple stellar populations with solar metallicity, changing the age 
and the dust content. The models are reddened according to a single parameter E{B — V) ~ that follows a Galactic 
reddening curve (e.g. lFitzpatricklll999| ). assuming R= Av/E{B — V) = 3.1 (other reddening laws will not introduce 
significant differences). The observed V-i and i-H colours constrain the stellar M/L in the observed H band, which 
is used to determine the stellar mass content, following the standard methodology of this paper. The panels, from 
top to bottom, show the best SSP-based luminosity-weighted age, stellar mass, and as a function of the reddening 
parameter E(B — V). We find that dust "conspires" with age such that an increase in dust is compensated by a 
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["] [kpc] 



Ml(< 2Rc) Ms(< 2R,o) 



["1 



nlO 



[kpc] 



Q0047 0.485 3. BO 2.20 12.82 0.880 ± 0.025 1.45 ±0.04 l.Sl^Jj Jj- 33.273| g'^ 11.58 ± 0.43 5.04^|lj-^j^ 189. 23^-^ 538.4 G(9) 

Q0142 0.490 2.72 2.23 13.10 0.703 ± 0.013 2.64 ± 0.05 2.54+^ T^J 45.81 + ^-^i, 13.45 ± 2.05 2.92 + iJ-^^ 245.8 + ^^-^ 666.5 

— 0.83 —5.50 —0.69 —46.8 

MG0414 0.960 2.64 2.08 16.01 0.734 ± 0.018 1.85 ± 0.05 1.67+^'^^ \02.7Q'^\^^. 19.90 ± 2.29 5.62 + -f-^^ 247. O^'^^^^^ 635.0 1 



B0712 0.410 1.34 1.29 6.82 0.702 ±0.016 1.15 ± 0.03 1.02 + ° °^ 1^1312 04 5.46 ± 0.52 3.27+JJ-^2 164.1 + g-^| 295.8 1 

HS0818 0.390 3.21 2.56 13.15 0.679 ± 0.013 3.27 ± 0.06 2.51 + ^5'^^ 36.88 + ^ |^^ 9.44 ± 0.97 5.42 + J-Jj^ 245.8 + ^*^1 517.8 1 

RX.J0911 0.769 2.8 3.22 23.16 0.725 ±0.011 3.09 ± 0.05 2.27+^-^|^ 73.04+^ 22 14.52 ± 1.87 2.99 + q"!^^ 27S.6 + g Q 590.4 C 

BRI0952 0.632 4.5 0.99 5.27 0.619 ± 0.021 1.04 ± 0.04 0.82 + °^^^ 14.76 + |-^^ 4.98 ± 0.35 2.78+55"^* ^"^^ ■'^t.^o% 233.5 G(5) 

Q0957 0.356 1.41 6.17 29.98 1.491 ± 0.018 3.51 ± 0.04 2.39+°-^li 151.29+^1:^ 20.92 ± 0.96 2.52+!!-?n ^''^'-slli'n 950.0 C 



LBQS1009 0.880 2.74 1.54 11.56 0.963 ± 0.028 1.27 ±0.04 l-97+^-3g 64.73 + 24 46 6.62 ± 0.96 2.14+JJ-g^ 220.9 + g5 g 252.1 

B1030 0.600 1.54 1.32 8.60 0.675 ± 0.019 2.05 ± 0.06 2.11 + °-|g 55.09 + 3 g2 5.16 ±0.80 1.12+JJ-^^ 256.8 + ^^^g 277.2 

HE1104 0.730 2.32 3.19 22.54 0.681 ± 0.010 3.08 ± 0.05 3.11 + 9 5? 72. 80 + ^?-^^ 14.43 ±1.62 S.SO+i.'^^ 302.9 + ^^9 604.5 

— U.^-J — O.ZZ — U.oU — 11, U 

PG1115 0.310 1.72 2.43 10.76 0.478 ± 0.009 2.94 ± 0.06 2.50+9 ?^ 16.80 + ^ 5.61 ± 1.19 3.68+i,-^?i^ 354.0 



B1152 0.439 1.02 1.56 8.60 0.691 ± 0.013 1.62 ± 0.03 1.43 + °-lii 30.43 + ^-^5! 8.57 ± 0.65 2.84+^^-5^ 216.9+5?^-^ 431.4 

— 0.25 —5.86 —0.50 —37.6 



B1422 0.337 3.62 1.29 6.02 0.241 ± 0-003 4.49 ± 0.06 3.56+^-^3 5.72 + ^ 2^ 2.91 ± 0.45 6.40 + ^ ^2 162.9+^j^g 231.3 G{17) 



SBS1520 0.710 1.86 1.57 10.97 0.947 ± 0.018 1.28 ± 0.02 0.95 + °°^ 41.98+j|2 11.49 ±1.41 2.63 + p^g 198.3+gj 434.3 G(4) 
B1600 0.420 1.59 1.38 7.45 1.015 ± 0.007 1.00 ±0.01 0.75 + ^ ^^ 16.38 + '^ g^ 3.93 ± 0.17 6.37+^-|^ 154.4+'^2^1 236.0 G(6) 



31608 0.630 1.39 2.10 13.92 0.839 ± 0.047 1.82 ±0.01 I '^S + g ^Q 42.49+2g'52 27.99 ± 1.63 2.44+q-^| 266.8+*^^"g 972.9 G(8) 



MG2016 1.010 3.3 3.36 26.22 0.406 ± 0.009 6.12 ± 0.14 ll.l + 2'g 52.05 + ^'*^3^ 6.34 ± 1.99 0.89+Q"2g 308.6 + 26^5 242.6 C(69) 

B2045 0.870 1.28 1.93 14.46 0.950 ±0.019 1.48 ± 0.03 1.23 + ° ]^^ 173.07+21°2 14.05 ± 1.03 2.47+jJ-^2 338.8 + ^5)'^ 517.2 



HE2149 0.603 2.03 1.71 10.02 0.531 ± 0.008 2.59 ±0.04 1-673q 23 26.82^2 19 4.28 ± 0.47 0.79_j^"^g 191. 2^^^ 242.6 

Q2237 0.039 1.7 1.83 1.40 1.090 ± 0.014 0.89 ± 0.01 0.81 + !! !!} 2.76 + n ?| 1.15 ± 0.12 4.39+J'™ 145.2 + ^ * 212.4 



TABLE 1 

Full set of gravitational lenses used for this analysis. All quantities in the table assume Hq = 72 km s~-^ Mpc"-', = 0.3 
AND = 0.7. The underlined values show maximum and minimum. A6 is the image separation. For systems with more than 
TWO images the maximal image separation between two images is given. Columns Re, Rleks/Re and REm/RE contain Petrosian 

radii determined in the observed ff-BAND WITH la ERROR BARS. NOTE THAT He™ IS COMPUTED FROM THE CRITICAL SURFACE DENSITY 
OF THE PIXELATED MAPS. THE TOTAL AND STELLAR MASSES ENCLOSED WITHIN 2Re ARE GIVEN IN THE FOLLOWING TWO COLUMNS. THE 
STELLAR MASS-TO-LIGHT RATIOS IN THE REST-FRAME V-BAND ARE MEDIAN VALUES OF ALL MODELS. CTeeks DENOTES THE VELOCITY 
DISPERSION AT Reeks- COLUMN i?vm GIVES THE VIRIAL RADIUS CALCULATED USING EQ. [2] AND OUR STELLAR MASS VALUES IN COMBINATION 
WITH THE Afs-TO-Afhalo RELATION FROM GUO FT AL. (2010). THE COLUMN LABELED Env CONTAINS ENVIRONMENTAL INFORMATION. "C" 
DENOTES A CLUSTER ENVIRONMENT, "G" A GROUP ENVIRONMENT AND "1" A LENS WITH ONLY ONE KNOWN COMPANION. If KNOWN THE 
NUMBER OF GROUP MEMBERS IS GIVEN IN PARENTHESES. REFERENCES FOR GIVEN VALUES ARE MENTIONED IN SECTION [2] COLOURS AND 
MAGNITUDES ARE IN AGREEMENT WITH COMPARABLE QUANTITIES IN RUSIN FT AL. (2003). 



younger age to give the same colours, yielding a small variation of the estimated stellar mass with respect to dust 
reddening. Most importantly, the value of worsens for high amounts of reddening. Hence we can safely say that 
even in the case of B1600 and B1608, the systematics on the stellar mass cannot be any larger than about 0.3 dex in 
log(Ms), shown as a shaded grey region in the middle panels of the figure. The other lens from our sample that could 
be affected by dust, Q2237 (i.e. the bulge of a late-type galaxy) is at a very low re dshift (z = 0.039), s o that stellar 
masses are determined from rest-frame H-band, which is even less sensitive to dust (jFerreras et al.|[2C)10|) . 
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TABLE 2 

List of lens systems and how they were treated to obtain surface bric;htness profiles of the lens galaxies. From left to 
RiciiiT the Lens-ID, the number of lensed images, the bands for which fitting was feasible, the PSF picking method, fitted 
and masked objects and necessary constraints are given. The column PSF includes the abbreviations tt for a PSF created 
WITH TinyTim, A for the outermost and thus fairly isolated quasar image, *(0414) for an isolated star close to lens 
MG0414, it for the iteration method applied to the most isolated image and A(1030) for an isolated quasar image taken 
from lens B1030 used for subtracting quasar images and for the convolution of the whole image. In column fitting WE 
summarize the number of objects, not significantly contributing to the lens mass, which are fitted WITH previously picked 
PSFs (P) AND Sersic profiles (5). In column masking xP refers to a number of x masked out point sources, mostly quasar 
images, 6 denotes resolved but indistinct objects which are not necessarily connected to the lens mass and henceforth 
excluded from each fit. Point sources like foreground stars indicated by * are also masked out. A,B refer to quasar 

images which could be masked out given their separation from the lens galaxy. The last column states the type of 
constraint used if necessary to prevent each fit from diverging. The constrained parameters ARE: effective radius Re, 
Sersic index n, axis ratio h/a, magnitude mag, position angle PA, sky background S and position of the Sersic profile xy. 

Parameters fixed at a certain value are underlined. 
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object B0712+472 
redshifts 0.41 1.34 
shear -45 

-0.013 -0.804 
0.795 -0.156 
0.747 -0.292 

-0.391 0.307 



object 1030+071 
redshifts 0.600 1.540 
shear +45 

0.846 1.097 
-0.085 -0.156 




object B1152+200 
redshifts 0.439 1.019 
shear +30 
double 

0.549 0.978 
-0.387 -0.268 



object B1422+231 
redshifts 0.337 3.620 
shear +45 

1.079 -0.095 
0.357 0.973 
0.742 0.656 
-0.205 -0.147 




object B1600+434 
redshifts 0.42 1.59 
shear +60 

0.610 0.814 
-0.110 -0.369 51 



object B1608+656 
redshifts 0.630 1.394 
shear -30 

-1.300 -0.800 

-0.560 1.160 31 

-1.310 0.700 5 

0.570 -0.080 40 



object B2045+265 
redshifts 0.870 1.280 
shear +30 

0.824 
0.035 
0.576 
.183 



1.121 
1.409 
1.255 
-0.507 ■ 




object BRI0952-012 
redshifts 0.632 4.500 
shear -45 
double 

-0.396 0.506 
0.300 -0.203 




object HE1104-181 
redshifts 0.730 2.320 
shear +30 
double 

-1.927 -0.822 
0.974 0.510 161 



% 



object HE2149-275 
redshifts 0.603 2.03 
shear +45 

0.736 -1.161 
-0.173 0.284 103 




object HS0818+123 
redshifts 0.390 3.115 
ptmass -2 2 2 10 

-1.657 -1.475 
-0.152 0.592 




object LBQS1009-025 
redshifts 0.880 2.740 
shear +90 
double 

-0.537 1.097 
0.133 -0.286 



object MG0414+053 
redshifts 0.960 2.640 
shear -30 
- 3 quads - 




object MG2016+112 

zlens 1.01 

shear +45 

- 3 ptmasses - 

multi 3 3.27 

-1.735 1.778 1 
0.317 -1.470 1 
1.269 0.274 2 




object PG1115+080 
redshifts 0.311 1.722 
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-0.968 -0.105 146 
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Fig. 15. — PIXELENS input data and projected mass distribution for the lens sample. The black dots mark the multiple images. All 
maps have a radius of 15 pixel. All mass maps have a radius of 2 Ricns, which corresponds roughly to 2 Rsin- All lens properties as well 
as respective length specifications are in Table [T] 
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